Discovery of High-Latitude CO in a HI Supershell in NGC 5775 
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ABSTRACT 

5/3 , We report the discovery of very high latitude molecular gas in the edge-on spiral galaxy, 

NGC 5775. Emission from both the J=l-0 and 2-1 lines of 12 C0 is detected up to 4.8 kpc 
away from the mid-plane of the galaxy. NGC 5775 is known to host a number of HI supershells. 
The association of the molecular gas (Mh 2 ,F2 = 3.1 xlO 7 M Q ) reported here with one of the HI 
supershells (labeled F2) is clear, which suggests that molecular gas may have survived the process 
which originally formed the supershell. Alternatively, part of the gas could have been formed 
in situ at high latitude from shock-compression of pre-existing HI gas. The CO J=2-l/J=l-0 
line ratio of 0.34±40% is significantly lower than unity, which suggests that the gas is excited 
subthermally, with gas density a few x 10 2 cm~ 3 . The molecular gas is likely in the form of 
cloudlets which are confined by magnetic and cosmic rays pressure. The potential energy of the 
gas at high latitude is found to be 2 x 10 56 ergs and the total (HI + H2) kinetic energy is 9x 
10 53 ergs. Based on the energetics of the supershell, we suggest that most of the energy in the 
supershell is in the form of potential energy and that the supershell is on the verge of falling and 
returning the gas to the disk of the galaxy. 

Subject headings: galaxies: halos — galaxies: ISM — galaxies: individual (NGC 5775) — galaxies: spiral 
- ISM: bubbles 

1. Introduction 
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studied by Heiles (1979, 1984). The large energies 
found in the supershells (> 10 53 ergs) imply that 
these structures must have a tremendous influence 
on the structure of the interstellar medium. In ad- 
dition, supershells that break through the gaseous 
disk to reach high galactic latitudes may be a 
source of star formation in the halo. For example, 
supershells may act as "chimneys" through which 
hot gas from the disk funnels to the halo (e.g., the 
"chimney model", Norman & Ikeuchi 1989). This 
hot gas may cool and eventually form stars in the 
halo. Alternatively, the molecular gas in the su- 
pershell may reach high-latitude and directly pro- 
vide raw material for star formation in the halo. 
Although which of these two scenarios is at work 
cannot be distinguished easily, the study of molec- 
ular gas in supershells still provides an important 
clue to high-latitude star formation. Molecular 
gas at high latitude also presents an important as- 
pect in the understanding of the global evolution 
of the interstellar medium in spiral galaxies. In 
the Milky Way Galaxy, the study of supershells 
is hindered by difficulties with distance determi- 
nation and the resulting confusion. In external 
galaxies, these problems are minimized. 

NGC 5775 is an edge-on (i = 86°), infrared- 
bright (Lfir = 2.6xl0 10 L Q ) galaxy at a distance 
of 24.8 Mpc. Irwin (1994) observed this galaxy 
and its face-on neighbour NGC 5774 in HI using 
the VLA and provided models for their HI distri- 
butions. She showed that the two galaxies may 
be engaging in an early phase of interaction, with 
two HI bridges connecting them. Numerous HI 
arcs and extensions beyond the disk of NGC 5775 
are also observed. Six HI supershells were cata- 
loged by Lee (1998). In a multi- wavelength study 
of NGC 5775, Lee et al. (2001) report spatial 
correlations of HI, radio continuum, X-rays and 
far-infrared emission at the positions of the three 
largest HI supershells, labeled Fl through F3 in 
Fig. 1. 

2. Observations 

Observations were performed with the IRAM 
30m telescope 2 on 2000 October, 2001 April and 
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Right Ascension (J2000) 



Fig. 1. — HI column density contours superim- 
posed on the Digital Sky Survey image of the 
vicinity of NGC 5775. The white crosses indicate 
the centres of NGC 5775 (left centre), NGC 5774 
(upper right) and IC 1070 (bottom). Contour 
levels are at 1, 5, 10, 17.5, 25, 40, 60, and 
82.6xl0 20 cm~ 2 . The HI supershells are labeled 
Fl to F3. 
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November. The CO J=l-0 line at 115 GHz and 
J=2-l line at 230 GHz were observed simultane- 
ously using two dual mixers tuned in single side- 
band mode. The observations were done in "wob- 
bler" switching mode with a wobbler throw of 
240", resulting in very flat baselines. Relevant ob- 
serving and spectral parameters are listed in Ta- 
bic f. Data reduction was carried out using the 
GILDAS software package. Gaussian fits were ob- 
tained in order to find the peak antenna tempera- 
tures, central velocities, full-width-half-maximums 
(FWHMs) and the integrated intensities of the 
lines. For clarity of presentation, we define the 
x and y directions to be parallel and perpendic- 
ular to the galaxy's major axis, respectively. All 
offsets are given in the form (x-offset,y-offset) in 
units of arcseconds, where the offsets are measured 
with respect to the observed centre at a(J2000) 
= 14 h 53 m 57!7, <5(J2000) = 3°32'40'.'0. Note that 
x-offset>0 represents the eastward direction and 
y-offsct>0 represents the northward direction. 



Table 1: Observing & Spectra Parameters. 



Parameter 


115 GHz 


230 GHz 


T sys a (K) 


250 


350 


HPBW b (") 


21 


11 


B e// C 


0.75 


0.52 


F e// d 


0.95 


0.91 


Channel Width (kms -1 ) 


26 


21 


Spectra rms (mK) 


2 
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"Typical system temperature. 
b Half-power-beamwidth. 
c Main beam efficiency. 
d Forward efficiency. 



3. Results 

We initially searched for CO emission at a few 
selected positions in all three HI supershells (see 
Fig. 1). Although all three showed tentative detec- 
tions, the spectral lines at F2 are the most obvious. 
We therefore proceed to map the CO emission in 
F2. The observed 12 CO J=l-0 and J=2-l spectra 
are presented in Fig. 2, superimposed on the HI 
total intensity map. Emission from both CO tran- 
sitions is detected up to 40" (4.8 kpc) away from 
the mid-planc of the galaxy. Within the region 
(35,-25) to (60,-40), roughly the size of the CO 
J=l-0 beam, the average peak values of Tmb for 
CO J = 1-0 and 2-1 arc about 13 mK (3.2 sigma) 
and 13.5 mK (2.3 sigma), respectively. The fact 
that there is an absence of CO emission near the 
HI "hole" at the offset of (55,-30) (see Fig. 2) sug- 
gests that the molecular gas mimics the shell-like 
distribution of the HI feature F2. 

Emission from a flared or warped disk may re- 
semble gas at high-latitude when seen in projec- 
tion. However, in such a case the line-of-sight cen- 
tral velocity of the emission line, which presum- 
ably originates from the outer part of the disk, 
would be much closer to the systemic velocity of 

(IGN). 




Fig. 2. — (a) Observed CO J=l-0 spectra superim- 
posed on the HI column density map of F2. The 
HI contour between x-offsct = 50" and 55" and 
y-offsct = -35" and -30" represents a decrease in 
column density. Horizontal axes of the spectra are 
from -500 kms" 1 to 500 kms -1 and vertical axes 
are plotted from -7 mK to 40 mK in the Tmb 
temperature scale, (b) As in (a) but for CO J=2- 
1. Horizontal axes of the spectra are plotted from 
-300 kms -1 to 300 kms -1 and the vertical axes 
from -7 mK to 28 mK. 
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the galaxy. We therefore compare the average cen- 
tral velocity and average FWHM of all CO 1-0 
spectra in F2 to a spectrum along the major axis. 
The average central velocity and average FWHM 
at F2 are 131±19 kms" 1 and 83±23 kms" 1 , re- 
spectively. The uncertainties are the standard de- 
viations of all the spectra and are both smaller 
than the velocity resolution. At (40,0) (spec- 
trum not shown), the corresponding values are 
126 kms -1 and 69 kms -1 , respectively. The cen- 
tral velocities comparison shows that the gas in 
F2 shares the same circular velocity as the major 
axis gas at the same galactocentric distance and 
does not originate from the outer disk. The nar- 
row line- widths support this, since gas from the 
outer disk would have wider spectra. Therefore, 
we conclude that the gas associated with F2 is in- 
deed at high galactic latitude and not from a flared 
or warped disk. Detailed analysis of the kinematic 
structures of the galaxy and the supershell will be 
forthcoming. 

In order to study the molecular gas that truly 
belongs to the supershell at high latitude, we have 
to eliminate the contribution of CO from the disk. 
We assume that the vertical CO distribution of 
the galaxy can be represented by a single gaus- 
sian which peaks at the mid-plane, and that the 
emission from the supershell at high latitude is 
superimposed on the wing of the gaussian. The 
emission from the supershell can then be isolated 
by subtracting the gaussian model from the ob- 
served data. To fit the gaussian, we average the 
CO integrated intensity map from x-offsets = 30" 
to 65", to increase the signal-to-noise. A gaussian 
was fitted from y-offsct = +20" to -20" (hence 
docs not include the high- latitude emission). The 
results of the fits are good, the residuals being less 
than 5% of the peak values in both transitions. 
At the centre of the supershell (y-offset = -32.5") 
, the difference in integrated intensities between 
the data and the model gives the CO emission 
from the supershell alone. The total (disk con- 
tribution plus supershell contribution) integrated 
intensities within a 21" beam at the centre of F2 
are 1.28 and 1.05 K-kms" 1 for the 1-0 and 2-1 
transitions, respectively, where the temperatures 
are in the Tmb scale. The uncertainties in these 
values, taking into account noise, calibration and 
baseline errors, are about 30%. For proper com- 
parison, the CO J = 2-1 data were smoothed to the 



resolution of the 1-0 data. After subtracting the 
disk contribution, the CO 1-0 and 2-1 integrated 
intensities are given by Iio = 0.89 K-kms -1 and 
I21 = 0.30 K-kms -1 , respectively. The line ratio of 
the two transitions, I21/T07 is therefore 0.34±40%. 

For comparison, we obtained the line ratio at 
mid-plane corresponding to the same galactocen- 
tric radius of F2. We use our JCMT CO J = 2-1 
integrated intensity (= 17.4±20% K-kms -1 ) at the 
position equivalent to (40,0) (Lee 1998), and the 
IRAM 30m CO J = 1-0 integrated intensity (= 
20.2±20% K-kms- 1 ) at the same position. The 
JCMT beam has the same HPBW as the IRAM 
30m beam at these frequencies. The resultant line 
ratio at the mid-plane is 0.86±30%. 

The total molecular gas mass in F2 (exclud- 
ing the disk contribution), obtained using the CO 
J=l-0 integrated intensity, is found to be Mh 2 ,f2 
= 3.1 xlO 7 M Q within the 21" beam (= 2.5 kpc in 
diameter) , assuming a CO-H2 conversion factor of 
3xl0 20 cm^K-kms -1 (Young & Scoville 1991). 
Since the CO and the HI distributions within F2 
agree fairly well, we have probably detected most 
of the CO emission in the supershell (see Fig. 2). 

We have also obtained the HI integrated inten- 
sity within a 21 "-beam at F2 following the same 
procedure outlined above. The HI data have been 
published in Irwin (1994) and in Lee et al. (2001). 
The HI integrated intensity, after subtracting the 
disk contribution, is 0.11 Jy/beam-kms -1 and the 
corresponding HI mass is Mhi,f2 = 1.6xl0 7 M Q , 
or 50% of Mij 2)F2 . The total gas mass (HI + H 2 ) 
in F2 is therefore 4.7xl0 7 M Q . 

4. Discussion 

This paper reports the highest latitude (4.8 kpc) 
molecular gas detected to date in any galaxy. In 
our own Galaxy, high-latitude molecular clouds 
are found in association with the HI shells around 
two OB associations, the Per OB3 and the Sco 
OB2 associations (Bhatt 2000); and in association 
with a hot bubble between Ccpheus and Cas- 
siopeia (Grenier et al. 1989). There are only a 
few examples of detections of discrete extraplanar 
molecular gas in external galaxies. In NGC 891, 
a molecular spur extending to 0.5 kpc above the 
mid-plane was detected by Handa et al. (1992). 
Possible formation mechanisms considered for this 
structure include ejection by a superbubble or the 
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buoyancy force on the gas due to Parker instabil- 
ity. Recently, molecular gas up to 3 kpc above the 
mid-plane was discovered in the molecular out- 
flow in M 82 (CO detection by Taylor, Walter, & 
Yun (2001); SiO detection by Garcfa-Burillo et al. 
(2001)). Many of these examples of discrete high- 
latitude molecular gas are thought to be related 
to star formation activities in the disk such as 
stellar wind or supernovae swept-up shells. At the 
highest latitudes, the molecular gas may represent 
the component of the shells that broke through 
the thin disk to reach the halo. 

4.1. The Physical Conditions in the Su- 
pershell 

The line ratio (0.34±40%) for the extraplanar 
gas obtained is significantly lower than the corre- 
sponding ratio at the mid-plane. Even taking into 
account the uncertainty, its highest possible value 
is less than 0.5. Assuming that the gas is thermal- 
ized so that the excitation temperature of all rota- 
tional levels is the same, the interpretation of such 
a low line ratio is that the excitation temperature 
is low (T ex ~ 5 K). Although no estimate of the 
gas or dust temperature in the supershell is cur- 
rently known, the possibility of shocks, deduced 
from the trend of radio continuum spectral indices 
(see Lee et al. 2001), suggests that the gas kinetic 
temperature may be > 5 K due to shock heating. 
In this case, the gas density must be low (~ a few 
hundred particles per cm 3 ) so that thermalization 
does not occur. Using a Large Velocity Gradient 
(LVG) analysis, we find that such a low line ratio 
is consistent with a gas density n# 2 < 1000 cm~ 3 , 
although wc cannot constrain the gas kinetic tem- 
perature. It is noted that the line ratio is expected 
to be closer to one in shock regions due to the 
elevated gas kinetic temperature, which tends to 
populate the upper rotational levels. Better con- 
straints of the physical parameters will have to 
await data of higher CO transitions. 

Given the gas density and gas mass within the 
21" beam, the filling factor of the molecular gas 
is of order 10~ 3 . That is, the molecular gas must 
be in the form of small cloudlets. We can esti- 
mate whether these cloudlets are in pressure equi- 
librium with the surrounding environment. The 
ambient pressure surrounding the cloudlets comes 
from the magnetic pressure (P m ), cosmic rays 
pressure (Pen), and the gas pressure. The first 



two are equal in the case of energy equipartition 
and are both given by B 2 /87r, where B ~ 3 /uG 
is the equipartition magnetic field (Duric, Irwin 
& Blocmcn 1998). The gas pressure is the com- 
bined ionized (Pi) and neutral gas pressure (Phi), 
each given by P = nkT, where n is the gas num- 
ber density, k is the Boltzmann constant and T 
is the kinetic temperature of the gas. We assume 
the ionized gas is at T = 10 4 K and the electron 
density n = 9xl0~ 4 cm~ 3 at 4 kpc above mid- 
plane (from the electron density distribution given 
in (Collins et al. 2000)). For the HI gas we use T = 
8000 K and density n HI = 2xl0~ 3 cm~ 3 at 4 kpc 
above midplane (from the vertical density distri- 
bution given in Irwin (1994)). From these values, 
we obtain P m ~ Pen. ~ 4xl0~ 13 erg cm -3 , Pj = 
2xl0~ 15 erg cm~ 3 , and P H i = 2xl0~ 15 erg cm~ 3 . 
The averaged internal pressure of the molecular 
gas is P doud = lxl0~ 12 erg cm -3 , using parame- 
ters given in the previous paragraph. It seems that 
the magnetic and cosmic rays pressure combined 
is sufficient to confine the clouds without internal 
gravitation. 

4.2. Origin of Molecular Gas at High Lat- 
itude 

Observations of molecular gas at high-latitude 
are interesting in that they may explain high- 
latitude star formation. The key question is: how 
does the molecular gas reach such high latitude? 
Two possibilities may be considered. First, the 
molecular gas was formed in the disk and was 
ejected by the outflow that formed the HI super- 
shell; and second, the molecular gas was formed at 
high-latitude. Of course, a combination of these is 
also possible. 

The gravitational potential energy (Epe) of a 
cloud at a height, z, above the galactic plane is 
calculated by solving the Poisson equation for a 
disk flattened in one-dimension. We have assumed 
the galactic value for the stellar scale height and 
the solar neighbourhood value for the stellar mass 
density at mid-plane (325 pc and 0.175 M pc -3 , 
respectively, from Freeman (1987)) since the cor- 
responding values for NGC 5775 are not known. 
Including both HI and H 2 gas, E PE = 2x 10 56 ergs 
in F2. Note that Pipe scales as the square of the 
stellar scale height so that if the scale height of 
NGC 5775 is twice that of the Milky Way, P PE 
will increase by a factor of 4. We estimate the to- 
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tal kinetic energy (Eke) of the gas (HI + H2) 
to be Eke = 9xl0 53 ergs. This includes the 
kinetic energy from the HI gas in F2, which is 
found to be 6xl0 53 ergs, using Mhi,f2 calcu- 
lated above and assuming the expansion veloc- 
ity is 62.5 kms" 1 (Lee et al. 2001); and the ki- 
netic energy from the H2 gas, which is found to 
be 4xl0 53 ergs. The expansion velocity of the 
H 2 gas is taken to be half the averaged FWHM 
of all the CO J=l-0 spectra (see § 3) in F2 mi- 
nus, in quadrature, the expected linewidth from 
differential rotation within the beam (44 kms -1 ) 
and the channel width (26 kms -1 ), and is equal 
to 32.5 kms -1 . The predicted differential rota- 
tion within the beam is calculated using the HI 
rotation curve model given in Irwin (1994). As- 
suming that the supershell expands isotropically 
(a fair assumption based on the near circular mor- 
phology of F2), then the line-of-sight expansion 
velocity observed is a good estimate of the true 
expansion velocity of the supershell. Therefore, 
Eke of the supershell is more than two orders of 
magnitude lower than Epe- The large difference 
in Epe and Eke suggests that the supershell is 
being observed at a time when it is about to fall 
back to the disk so that most of the energy is in 
the form of potential energy We can therefore vi- 
sualize an explosion at the mid-plane that created 
an expanding bubble of hot gas. In the process, 
neutral gas (atomic and molecular) is entrained in 
the ejecta, forming the supershell, reaching large 
z-height. Under the influence of the gravitational 
potential, the supershell is now about to plunge 
back, returning the gas to the star-forming disk. 
Such a recycling process has been proposed in the 
galactic fountain model (Bregman 1980). 

Another possibility for the molecular gas to ex- 
ist at high latitude is that H 2 gas can be formed 
in situ via shock-compression (e.g., Magnani 1987; 
Elmegreen 1988), occurring within HI pre-existing 
in this region. The effect of a strong shock prop- 
agating through the interstellar medium is that 
dense molecular cloudlets can form due to thermal 
instability in the shock regions (Koyama & Inut- 
suka 2000). The existence of Ha emission (possi- 
bly due to shock-ionization, see Lee et al. 2001) in 
F2 suggests that at least some of the molecular gas 
may be formed this way. It is more difficult to ex- 
plain how the HI come to reside at high-latitude, 
in view of the large mechanical energy required 



(equals to Epe)- Various mechanisms have been 
examined in Lee et al. (2001), including supernova 
explosions and cloud impacts. In any case, this 
material would likely have originated within NGC 
5775 since the mean radial velocity agrees with the 
orbital velocity of the disk. Finally, assuming typ- 
ical Galactic dust-to-gas ratio, the high-latitude 
molecular clouds could be emitting substantial in- 
frared radiation and has indeed been observed at 
850yum (Brar et al., private communication). 

5. Conclusions 

This paper reports the detection of high- 
latitude molecular gas in NGC 5775. The shell-like 
distribution of the CO emission coincides exactly 
with that of the HI supershells, suggesting that we 
have detected the molecular shell associated with 
the HI supershell. The existence of the molecular 
shell means that molecules are not destroyed dur- 
ing the ejection of the supershell and is entrained 
in the expanding flow to reach high latitude. Some 
of the molecular gas may have been formed in 
situ, via shock-compression of pre-existing HI gas. 
The CO J=2-l/l-0 line ratio (0.34±40%) suggests 
that the gas density in the supershell is low and 
the gas is subthermally excited. The molecular 
gas is probably in the form of cloudlets which are 
confined by magnetic and cosmic rays pressure. 
Based on energetics ground, we proposed that the 
supershell may be at a stage where it is about to 
plunge towards the disk of the galaxy, returning 
the gas to the bulk of the gas reservoir of the 
galaxy. 
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